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摘 要 : 选择 草地 类 型 丰富 多 样 的 北 疆 作 为 研究 区 。 基 于 MODIS NDVI 数 据 , 获 取 时 间 累 积 归 一 化 植被 指数 TI- 
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NDVI 和 年 最 大 NDVI(NDVImax) , {ff BY GIS 空间 分 析 技 术 , 25 


异 系数 (Cv) 及 Mann-Kendall 非 参数 趋势 检验 等 方法 ,对 


2000 一 2019 年 北 疆 地 区 草地 覆盖 动态 变化 进行 了 分 析 , 并 探究 了 TI-NDVI 和 NDVImax 对 草地 时 空 异 质 性 的 表达 能 
力 的 比较 优势 。 结 果 表 明 :(1) 用 NDVImax 和 TI-NDVI 表 征 的 北 疆 草地 呈现 明显 海拔 分 异 。TI-NDVI 总 体 随 NDVI- 


max 的 增 大 而 增 大 ,但 NDVImax 或 TNDVI 相 同 的 区 域 ,其 TI-NDVI 或 NDVImax 却 存在 较 大 差异 ,(2) 2000—2019 


年 北 疆 地 区 草地 TI-NDVI 和 NDVImax 总 体 呈 显著 增加 趋势 (P<0.01) ,但 草地 TI-NDVI 变 化 的 空间 分 异 与 NDVImax 
明显 不 同 , 全 区 17.55% 的 草地 TI-NDVI 变 化 趋势 与 NDVImax 变化 趋势 相反 。 尤 其 阿尔 奈 山 与 伊犁 河谷 ,高 覆盖 草 
地 分 布 区 的 NDVImax 与 TI-NDVI 均 呈 相 反 变 化 趋势 ,(3) 在 时 间 和 空间 维度 上 , 北 疆 山 区 高 覆盖 草地 TI-NDVI 的 C， 
均 高 于 NDVImax。TLNDVI 对 高 覆盖 草地 的 时 空 异 质 性 反映 更 敏感 ,能 在 一 定 程度 上 削弱 草地 动态 评价 中 NDVI 光 


饱和 缺陷 的 影响 。 
关键 词 : 时 空 异 质 性 ; TI-NDVI; NDVImax; 草地 ; JU38 


气候 变化 和 人 类 活动 带 来 了 全 球 或 区 域 生态 
环境 的 改变 。 草 地 是 陆地 生态 系统 的 重要 组 成 部 
分 ,草地 变化 构成 了 全 球 变化 的 重要 组 成 部 分 ”。 
草地 是 人 类 生存 的 重要 自然 资源 和 环境 基础 ,通过 
对 草地 变化 研究 能 为 人 类 有 效应 对 全 球 变化 市 来 
的 环境 问题 提供 重要 参考 和 支持 。 而 准确 且 敏 感 
的 参考 指标 是 草地 变化 研究 的 基础 ”。 

异 质 性 是 地 表 生 态 系 统 的 基本 特征 ,也 是 评价 
地 表 徐 被 动态 要 回答 的 基本 问题 。 为 充分 理解 草 
地 动态 的 特征 或 规律 ,学 者 们 在 时 间 和 空间 维度 上 
揭示 了 关于 草地 异 质 性 的 更 多 细节 。 传 统 的 实 
地 调查 的 草地 监测 方法 由 于 耗 时 和 劳动 密集 型 ,很 
难 在 大 尺度 上 实现 对 草地 的 时 空 监测 “”。 由 于 遥 
感 技术 在 获取 数据 的 时 效 性 空间 性 及 成 本 等 方面 
具有 无 法 比拟 的 优势 ,是 目前 区 域 植 被 动态 监测 的 
主要 方法 及 数据 来 源 * 。 借 助 于 遥感 数据 计算 得 
到 的 植被 指数 是 量化 地 表 植 被 状况 的 有 效 和 实用 


收 稿 日 期 : 2021-11-20; 修订 日 期 : 2022-02-28 


方法 。 各 植被 指数 中 归 一 化 植被 指数 (NDVD) 与 光 
合 有 效 辐 射 .生物 量 、 植 被 生产 力 等 紧密 相关 ,是 目 
前 表征 植被 生长 状态 及 植被 覆盖 度 中 应 用 最 多 的 
植被 指数 '"。 日 多 家 机 构 或 研究 学 者 已 制备 了 NO- 
AA-AVHRR NDVI, SPOT- VGT NDVI, TERRA- MO- 
DIS NDVI &l GIMMS NDVI 等 多 种 覆盖 全 球 的 NDVI 
成 品 数 据 , 以 方便 使 用 。MODIS 数据 因 其 较 高 的 空 
间 分 辩 率 ,并 且 改 进 了 在 红 光 波段 和 近 红 外 波段 的 
波幅 及 辐射 定 标 等 技术 ,使 其 更 详尽 地 反映 地 表 植 
被 的 空间 差异 ,已 被 广泛 应 用 于 区 域 植被 的 动态 监 
myo, 

在 对 NDVI 时 间 序 列 数据 预 处 理 过 程 中 ,年 最 
大 值 合成 法 (MVC) 常 被 用 于 计算 年 内 NDVI 最 大 值 
(NDVImax) ,来 表征 单个 生长 季 内 植被 生长 达到 最 
好 时 的 状况 ,并 作为 植被 在 年 时 间 尺 度 上 的 状态 指 
标 拉 呈 。 然 而 ,NDVI 在 浓密 的 植被 冠 层 中 容易 达到 
饱和 ,敏感 性 降低 ,削弱 了 NDVI 对 植被 时 空 异 质 性 
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的 表达 能 力 。 同 时 植被 的 动态 变化 不 仅 存 在 于 植 
被 能 达到 最 佳 时 的 状况 上 ,也 体现 在 植被 生长 过 程 
的 年 际 差异 上 ,如 植被 生长 季 长 度 的 变化 等 。 时 间 
BFR NDVI (Time Integrated NDVI, TI-NDVI) BE E X. 
为 年 内 NDVI 变 化 曲线 与 时 间 轴 围 成 的 面积 或 一 定 
时 期 内 NDVI 的 累积 值 2 。 基 于 地 面 监测 和 卫星 
遥感 的 多 项 研究 表明 ,在 农作物 .灌木 .草地 和 和 森林 
等 多 种 生态 系统 中 TI-NDVI 与 植物 生产 力 具 有 很 强 
的 相关 性 ”小 。 基 于 此 ,TI-NDVI 常 被 作为 生产 力 
的 替代 指标 广泛 用 于 植被 动态 .生态 功能 及 环境 脆 
弱 性 等 研究 之 中 ?2 。 相 对 于 NDVImax,TI-NDVI 
不 仅 积 累 了 植被 的 NDVImax, 还 积累 了 植被 在 返青 
期 和 枯萎 期 的 NDVI, 揭 示 更 多 关于 植被 异 质 性 的 
信息 。 

鉴于 此 ,本 研究 选择 草地 分 布 广泛 且 类 型 丰富 
的 新 疆 北部 ( 北 疆 ) 作 为 研究 区 ,在 16 d 合 成 的 MO- 
DIS NDVI 时 间 序 列 数据 的 支持 下 ,借助 于 GIS 空 间 
分 析 技 术 以 及 Mann-Kendall 非 参数 趋势 检验 和 变异 
系数 (Cv) 等 数理 统计 方法 ,研究 北 疆 草地 时 空 动 态 ， 
分 析 NDVImax 和 TI-NDVI 在 表达 草地 异 质 性 方面 
的 优势 比较 ,以 期 为 北 又 地 区 乃至 全 球 草地 有 效 评 
估 和 管理 提供 参考 。 


1 研究 区 概况 


研究 区 域 位 于 新 疆 北部 ( 北 疆 ,42°8'24"~48°3'N、 
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80°36'36"~95°55/12"E) , JEJE FS PERAK, FAR RE SP 
于 -157~7077 m, 主要 以 山地 和 盆地 为 主 ( 图 1) ,年 
均 降水 量 200 mm 以 下 ,年 均 气 温 -4~9 "C ,属于 温带 
大 陆 性 干旱 半 和 干旱 气候 。 区 域 中 部 的 准噶尔 盆地 
被 天 山 .阿尔 泰山 以 及 准 嘲 尔 西部 山地 环 统 , 西 风 
气流 从 西南 部 的 伊犁 河谷 地 带 进 入 盆地 ,后 被 东部 
天 山 阻 隔 ,无 法 继续 东 进 ,水汽 停留 在 原 地 ,给 伊犁 
河谷 带 来 丰富 的 降水 。 北 疆 地 区 降水 量 空间 分 布 
极 不 均匀 ,山区 降水 量 远 高 于 平原 和 盆地 ,西南 部 
多 于 东北 部 。 伊 犁 河谷 天山 和 阿尔 泰山 年 平均 降 
水 量 超过 200 mm, 属 于 半 干 旱 区 ,其 他 区 域 属于 干 
早 区 。 人 研 究 区 域内 除了 额 尔 齐 斯 河 为 外 流 河 ,其 他 
河流 均 为 内 流 河 , 补 给 水 源 主要 为 大 气 降水 冰雪 
融 水 和 地 下 水 。 


2 材料 与 方法 


2.1 数据 来 源 与 预 处 理 

研究 用 到 的 NDVI 数 据 为 美国 NASA EOS 数据 
中 心 (https://ladsweb.modaps.eosdis.nasa.gov) 发 布 的 
MODIS MOD13Q1 产品 。MOD13Q1 数据 为 16 d 合 
成 的 时 间 序 列 数据 ,空间 分 辨 紊 为 250 m, 每 年 23 
期 ,时 间 序 列 是 2000 一 2019 年 。 在 数据 预 处 理 过 程 
中 ,对 NDVI 时 间 序 列 数据 进行 了 镶 骨 、 转 投影 、 重 
采样 .研究 区 裁剪 以 及 Savitzky-Golay JE YZ AR HE Ț 
用 于 年 际 变 化 分 析 的 NDVImax 数据 是 通过 对 每 年 
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图 1 研究 区 位 置 及 地 形 
Fig. 1 Location and topography ofthe study area 
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23 期 数据 的 最 大 值 合成 (MVC ) 处 理 获 得 , 即 用 年 内 
植被 生长 达到 最 好 状况 时 的 NDVI 作 为 年 尺度 上 的 
NDVI. 

dE SEE 5 8 319 321 B] at BA X Senti- 
nel-2 假 彩色 合成 影像 的 目 视 解 译 获 得 。Sentinel-2 
影像 选取 2018 年 6 一 9 月 植被 覆盖 达到 最 大 , 且 云 
E < 10% 的 高 质量 影像 。 假 彩色 合成 影像 像 元 大 小 
为 10 mm。 对 解 译 获 得 的 草地 分 布 矢量 数据 进行 了 
顶 格 化 处 理 ,作为 裁剪 NDVI 数 据 的 掩 膜 。 草 地 分 
布 数据 的 栅 格 化 处 理 及 NDVI 数 据 重 采样 处 理 过 程 
中 , 像 元 大 小 均 设 置 为 230 mx250 m, 以 确保 其 与 
NDVI 图 像 的 空间 匹配 。 同 时 NDVImax < 0.1 的 像 元 
被 作为 无 植被 区 进行 了 剔除 。 
2.2 计算 TINDVI 

TI-NDVI 通 常 被 定义 为 生长 季 内 NDVI 变 化 曲 
线 与 时 间 轴 围 成 的 面积 或 NDVI 的 累积 值 。 累 积 值 
计算 方便 日 简单 ,因此 TI-NDVI 利 用 生长 季 内 NDVI 
的 累积 值 表 示 , 即 : 


305 
TI- NDVI= >》NDVI，i= 65,81,97,…,305 
i=65 


式 中 :i 代表 NDVI 数 据 的 日 期 。MOD13Q1 产品 每 年 
23 期 数据 ,从 每 年 的 第 1d 到 第 353 d 每 16d 一 期 数 
Jg ,并 在 计算 过 程 中 对 NDVI<0 数 据 进 行 了 剔除 ， 
避免 植被 尚未 变 绿 或 已 经 变 黄 时 的 NDVIS5。 由 于 
北 疆 植被 生长 季 通 常 开始 于 3 月 初 ,结束 于 11 月 
初 ,为 了 进一步 降低 非 植被 的 噪声 ,消除 一 年 中 早 
于 65 d 、 晚 于 305 d 的 NDVI。TI-NDVI 计 算 取 3 月 初 
(第 65 d) 到 11 月 初 (第 305 d) 时 段 内 共 16 期 NDVI 
数据 。 由 于 NDVI 的 取 值 范围 为 0~1.0, 因 此 ,理论 
上 TLNDVI 的 取 值 范 围 为 [0,16]。 
2. NDVImax 与 TINDVI 的 比较 

分 别 对 2000 一 2019 4F. NDV Imax 和 TI-NDVI 时 
间 序 列 影像 求 平均 值 , 获 得 NDVImax 和 TI-NDVI 多 
年 平均 影像 ,对 NDVImax 和 TI-NDVI 所 揭示 的 草地 
空间 格局 进行 差异 性 比较 。 并 提取 NDVImax 多 年 
平均 影像 所 有 像素 值 以 及 其 所 对 应 的 TLNDVI 值 ， 
以 此 计算 NDVImax 与 TLNDVI 的 相关 系数 ,评估 两 
者 相关 性 。 

利用 Mann- Kendall 非 参数 趋势 检验 方法 ,对 
NDVImax 和 TI-NDVI 时 间 序 列 影像 进行 逐 像 元 计 
算 , 获 得 2000 一 2019 年 NDVImax 和 TI-NDVI 变 化 趋 
势 空间 分 布 图 ,分 析 NDVImax 和 TLNDVI 变 化 趋势 
的 空间 分 布 。 同 时 ,对 NDVImax 变化 趋势 图 与 TI- 


NDVI 变 化 趋势 图 进行 琶 置 操作 ,制作 两 者 变化 趋 
势 差异 的 空间 分 布 图 ,分 析 NDVImax 与 TLNDVI 变 
化 趋势 的 空间 差异 。 

变异 系数 (Cv) 用 于 衡量 数据 系列 中 数据 点 的 离 
散 程度 。Cv 越 大 ,说 明 数 据点 相对 平均 值 离散 程度 
越 大 。 在 用 NDVImax 或 TI-NDVI 等 指标 衡量 植被 
动态 时 ,所 选 指标 的 Cv 越 大 ,能 揭示 的 植被 波动 就 
越 大 ,能 揭示 植被 波动 的 细节 就 更 多 ,对 植被 动态 
也 更 加 敏感 。 因 此 ,分 别 在 时 间 和 空间 维度 上 对 
NDVImax 和 TINDVI 的 Cv 进 行 了 计算 。 时 间 维 度 
上 ,分 别 对 NDVImax 和 TI-NDVI 时 间 序 列 影像 进行 
逐 像 元 计算 ,获得 时 间 Cv(CCvw) 的 空间 分 布 图 。 空 间 
维度 上 ,分 别 在 NDVImax 和 TI-NDVI 多 年 平均 值 影 
像 数 据 的 基础 上 ,利用 15x15 像 元 的 滑动 窗口 ( 即 移 
动 窗 口中 包含 225 个 像 元 ) 进 行 逐 像 元 计算 ,获得 空 
lis] Cy(Cvs) 的 空间 分 布 图 。 利 用 Cv 和 Cvs 空间 分 布 
图 ,在 时 间 和 空间 维度 上 ,对 比分 析 NDVImax 和 TI- 
NDVI 对 草地 异 质 性 表达 的 敏感 性 及 两 者 的 比较 优 
势 。 在 计算 Cws 时 ,笔者 也 尝试 了 5x5 像 元 和 9x9 像 
元 的 移动 窗口 ,其 Cw 空间 分 异 与 1Sx15 像 元 移动 窗 
口 的 结果 基本 一 致 。 


3 结果 与 分 析 


3.4 北 疆 地 区 NDVImax 和 TILNDVI 空 间 格 局 

在 像 元 尺度 上 计算 2000—2019 4E 4b 5 NDVI- 
max 和 TI-NDVI 的 年 平均 值 ,表征 NDVImax 和 TI 
NDVI 的 空间 格局 (图 2)。 

从 图 2 可 以 看 出 , 北 又 地 区 NDVImax 和 TI-ND- 
VI 的 空间 格局 基本 相同 , 均 表 现 出 明显 的 海拔 分 异 。 
北 疆 地 处 欧 亚 大 陆 干 旱 半 干 旱 区 腹地 ,水 是 制约 植 
被 生长 的 关键 因素 ,高 程控 制 着 水 的 空间 分 布 *™。 
阿尔 泰山 山区 、 天 山 山脉 和 塔 城 兮 地 周边 山区 降水 
相对 丰沛 ,草本 植物 生长 良好 ,草地 NDVImax 和 TI- 
NDVI 分 别 >0.5 和 >5.0。 而 在 海拔 3500 m 以 上 的 山 
区 ,由 于 冬季 较 长 ,气温 较 低 ,植被 的 生长 受到 制约 ， 
草地 NDVImax 和 TI-NDVI 相 对 较 低 Y", 一 般 分 别 < 
0.4 和 <3.0 ,而 在 准噶尔 盆地 和 伊犁 河谷 的 低 海 拔 地 
区 ,由 于 降水 稀少 ,气候 炎热 ,干旱 制约 了 植被 的 生 
长 ,草地 NDVImax 和 TLNDVI 最 低 ,NDVImax<0.2， 
TINDVI<2.0。 

对 比 图 2a 和 图 2b 可 以 看 出 ,阿尔 泰山 北部 和 天 
山西 南部 的 高 山区 域 ,TI-NDVI 的 低 值 区 域 (TI-ND- 


202210.00018v1 


chinaXiv 


ChinaXiv 合 作 期 刊 
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多 年 平均 NDVImax 
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多 年 平均 TI-NDVI 
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图 2 北 疆 年 平均 NDVImax(a) 和 TI-NDVI(b) 的 空间 分 布 
Fig.2 Spatial distribution of yearly averaged NDVImax (a) and TI-NDVI (b) in northern Xinjiang 


VI<2.0) 明 显 多 于 NDVImax 的 低 值 区 域 (NDVImax< 
0.3)。 这 主要 是 因为 此 区 域 生长 季 短 ,即使 其 ND- 
VImax 与 其 他 山区 差异 较 小 ,其 累积 NDVI 却 相差 很 
大 。 表 明 TI-NDVI LU NDVImax 更 能 表达 高 山区 域 
草地 空间 分 布 的 异 质 性 。 
3.2 北 疆 地 区 NDVImax 和 TI-NDVI 年 际 变 化 特征 

图 3 显示 了 北 疆 整个 草地 NDVImax 和 TI-NDVI 
的 年 际 变 化 。Mann-Kendall 趋势 检验 表明 ,2000 一 
2019 年 北 疆 草 地 NDVImax 和 TI-NDVI 均 呈现 显著 
增加 趋 热 (P<0.01),Mann-Kendall 统 计量 Z 分别 为 
2.75 和 2.69。 年 际 波 动 上 , 除 2000—2001 年 、 
2005—2006 4E. .2014—2015 年 和 2018 一 2019 年 外 ， 
NDVImax 和 TILNDVI 的 波动 变化 也 基本 一 致 。 

在 空间 上 ,NDVImax 和 TILNDVI 的 变化 趋势 存 
在 明显 差异 ,但 总 体 上 均 表 现 出 明显 的 海拔 分 异 
(图 4a, 图 4b)。NDVImax 或 TLNDVI 呈 上 升 趋势 的 
区 域 主要 分 布 在 准噶尔 盆地 伊犁 河谷 和 塔 城 贫 地 


NDVImax 


图 3 2000 一 2019 年 北 疆 草 地 区 域 平均 NDVImax 和 TI- 
NDVI 的 年 际 变化 
Fig. 3 Inter-annual variations of the regional averaged 
NDVImax and TI-NDVI in the grassland of northern 
Xinjiang from 2000 to 2019 


的 洼地 。 下 降 趋 势 主要 分 布 在 天 山 山 脉 .阿尔 泰山 
山脉 和 塔 城 舍 地 周围 山区 。 经 计算 , 北 疆 20.87% 草 
原 的 NDVImax 是 下 降 趋 势 , 其 中 2.31% 达 到 显著 水 
平 (Z<-1.96)。NDVImax 呈 上 升 趋势 的 草地 面积 
79.13% , 达到 显著 水 平 (Z>1.96) 的 草地 面积 
19.24%。TILNDVI 呈 上 升 趋势 和 下 降 趋势 的 比例 分 
别 为 13.80% 和 86.20% ,达到 显著 水 平 (IZI>1.96) 的 
比例 分 别 为 0.48% 和 30.22%( 图 5a)。 

图 4c 是 NDVImax 和 TLNDVI 的 变化 趋势 的 县 
加 图 。 北 疆 大 面积 草地 NDVImax 和 TINDVI 呈 相 
反 变 化 趋势 (图 4c 中 ILD 和 D-T) , 占 比 达到 17.55% 
(图 $b) , 且 多 位 于 植被 覆盖 度 较 高 的 山区 。 尤 其 是 
阿尔 泰山 与 伊犁 河谷 周边 的 山区 ,草地 NDVImax 和 
TI-NDVI 变 化 趋势 的 差异 性 明显 不 同 。 阿 尔 泰山 山 
区 草地 变化 主要 表现 为 NDVImax 增加 而 TI-NDVI 
减 小 ,而 伊犁 河谷 周边 山区 则 主要 表现 为 NDVImax 
减 小 而 TI-NDVI 增 加 。 整 个 北 疆 地 区 ,2000 一 2019 
年 ,8.56% 的 区 域 NDVImax 和 TI-NDVI 同 时 下 降 (D- 
D) ,空间 上 主要 分 布 在 天 山南 部 的 部 分 地 区 ; 
73.89% 的 区 域 NDVImax 和 TI-NDVI 同 时 增加 (I-1)， 
主要 分 布 在 植被 覆盖 度 较 低 的 准噶尔 盆地 ;5.24% 
的 草地 NDVImax 增 加 而 TINDVI 减 少 (LD) ,空间 上 
主要 分 布 在 在 阿尔 泰山 脉 以 及 准噶尔 盆地 中 部 和 
伊犁 河谷 的 部 分 区 域 ;12.31% 的 草地 NDVImax 减 少 
而 TI-NDVI 增 加 (D-DD) ,空间 上 集中 分 布 在 伊犁 河谷 
和 塔 城 贫 地 周边 山区 。 
3.3 ”NDVImax 与 TI-NDVI 在 草地 时 空 异 质 性 表 
达 上 的 对 比 

时 空 异 质 性 是 评价 植被 动态 的 基本 问题 之 
一 。 通 过 计算 NDVImax 和 TI-NDVI 的 时 间 和 空间 
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注 :1.96 为 Mann-Kendall 检 验 P=0.05 显著 水 平 上 Z 的 临界 值 。D-I 
表示 NDVImax 减 少 ,TI-NDVI 增 加 ;I-D 表 示 NDVImax 增 加 ,TI- 
NDVI 减 少 ;D-D 表 示 NDVImax,TI-NDVI 均 减少 ;LI 表示 
NDVImax,TI-NDVI 均 增加 。 下 同 。 

图 4 NDVImax(a) 和 TI-NDVI(b)Mann-Kendall 非 参数 检 
验 统计 量 Z 值 的 空间 分 布 及 变化 趋势 天 加 (c) 

Fig.4 Spatial distribution of the statistic Z of the Mann- 
Kendall non-parametric test for NDVImax (a) and TI-NDVI 
(b), and the superimpose d trend image of NDVImax and 
TI-NDVI (c) 


Cv ,来 对 比 两 者 在 表达 草地 时 空 异 质 性 的 差异 中 ?1。 
图 6a 显示 了 NDVImax( 图 6al) 和 TINDVI( 图 
6a2) 的 Cvs 的 空间 分 布 及 其 差异 (图 6a3)。NDVI- 


60 F (a) = NDVImax 
55 um TI-NDVI 


百分比 /% 
cou Sas 


«-1.96 -1.96-0 0~1.96 >1.96 


max il TI-NDVI FY Cvs 的 空间 格局 几乎 相同 。Cvs 
(Cvs>0.3) 较 大 的 区 域 ,多 位 于 天 山 和 阿尔 泰 的 高 海 
d [X ,这 些 区 域 NDVImax 和 TI-NDVI 均 随 着 海拔 的 
升 高 而 快速 增 大 或 减 小 (图 2), 因 而 Cvs 均 较 大 。Cvs 
较 小 的 区 域 , 主 要 分 布 在 准噶尔 盆地 低地 或 天 山 、 
阿尔 泰山 半山 腰 , 这 些 区 域 NDVImax 和 TINDVI 在 
空间 上 的 变化 不 大 (图 2)。 然 而 ,54.18% 的 区 域 TI- 
NDVI 的 Cw 大 于 NDVImax 的 Cvs( 图 6a3), 空 间 上 主 
BEA} fi TE ND Vmax 通常 较 大 (NDVI>0.6) 或 沿海 拔 
梯度 变化 较 大 的 山区 ,以 及 小 部 分 植被 覆盖 度 最 低 
的 地 区 (NDVImax<0.2 fl TI-NDVI«1.0) (图 6a3 , 图 
2) ,这 说 明 TI-NDVI 更 能 表达 山区 植被 覆盖 度 较 高 
地 区 的 草地 空间 异 质 性 。 相 反 ,在 准噶尔 盆地 地 区 
的 低地 和 伊犁 河谷 地 区 , NDVImax 和 TI-NDVI 都 很 
小 (NDVI<0.4 和 TLNDVI<4.0) ,TLNDVI 的 Cus 多 小 
于 NDVImax 的 Cys (Al 6a) ,这 意味 着 NDVImax 在 表 
达 低 覆盖 草地 的 空间 异 质 性 方面 更 敏感 。 

图 6b 显示 了 2000—2019 年 期 间 NDVImax (图 
6b1) ffl TI-NDVI ( K] 6b2) Ct 的 空间 分 布 及 其 差异 
(图 6b3)。 相 较 于 NDVImax 和 TI-NDVI 的 Cvs 空间 
格局 ,NDVImax 和 TI-NDVI 的 Cw 的 空间 格局 差异 很 
大 (图 6b1, 图 6b2)。 对 于 NDVImax 而 言 ,在 准噶尔 
盆地 . 塔 城 盆 地 和 伊犁 河谷 的 洼地 的 Cx 明显 大 于 
山区 。 相 反 , 除 准 噶 尔 盆地 植被 覆盖 度 最 低 的 部 分 
地 区 外 ,低地 和 山区 TLNDVI 的 Cv 差异 不 大 , 均 处 
于 较 低 水 平 (0.05<Cvw<0.1)( 图 6b2)。 低 地 的 植被 覆 
盖 率 通常 比 山区 低 的 多 。 植 被 覆盖 度 高 的 地 区 ND- 
VI 容 易 饱 和 ,导致 高 覆盖 度 草地 NDVI 的 异 质 性 被 
低估 ,如 天 山 和 阿尔 泰山 区 的 草地 。 图 6b3 为 TI- 
NDVI 的 Cw 减 去 NDVImax 的 Cw 的 差 值 图 ,从 图 上 
可 以 看 出 ,植被 覆盖 度 较 高 的 山区 TI-NDVI 的 Cvt & 
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: | | 
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图 5 Mann-Kendall 检 验 (a) 和 羡 加 趋势 (b) 的 Z 范 围 比例 
Fig. 5 Proportions of Z ranges of Mann-Kendall test (a) and superimposed trends (b) 
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图 6 NDVImax 和 TI-NDVI 时 空 Cv 的 空间 分 布 及 其 差异 
Fig.6 Spatial distribution of the spatial and temporal Cy of both NDVImax and TI-NDVI, and their difference 


KF NDVImax 的 Cvt。 这 意味 着 TI-NDVI 可 以 弥补 
NDVImax 饱 和 度 的 不 足 , 更 能 表达 高 覆 善 草地 的 时 
间 异 质 性 。 然 而 ,低地 区 域 TLNDVI 的 Cw 比 NDVTL 
max [fJ Cys 要 小 。 


4 讨论 


4.1 NDVImax 与 TI-NDVI 的 关系 

NDVImax 是 植被 在 一 个 生长 季 所 能 达到 的 最 
大 NDVI, 而 TINDVI 定 义 了 一 个 生长 季 植 被 的 累积 
NDVI, 这 是 NDVImax 和 TLNDVI 的 根本 区 别 。 然 
ifi , 1588 NDVImax 和 TI-NDVI 的 空间 变化 基本 一 致 
(图 四 ,图 7a 也 显示 TI-NDVI 随 NDVImax 的 增加 而 增 
加 ,这 表明 NDVImax 高 的 地 区 TI-NDVI 总体 也 较 
Ko NDVImax 与 TI-NDVI 的 Pearson 相关 系数 达到 
0.89(P<0.01)。 因 此 ,NDVImax 与 TI-NDVI 在 表示 
植被 空间 分 布 方面 有 很 多 共同 之 处 。 


从 图 7a 中 可 以 看 出 , 随 着 NDVImax 从 0.1 变 化 
到 0.7,TI-NDVI 离 散 分 布 在 2 个 明显 不 同 的 范围 内 ， 
这 表明 具有 相同 NDVImax 的 区 域 通常 可 以 分 为 两 
部 分 , 即 TLNDVI 较 小 的 部 分 和 TLNDVI 较 大 的 部 
分 (图 7b1)。 同 样 ,TLNDVI 相 同 的 区 域 也 可 以 用 同 
样 的 方式 分 成 两 部 分 (图 7b2)。 表 明 在 相同 NDVI- 
max 或 TI-NDVI 的 区 域 可 能 有 完全 不 同 的 TI-NDVI 
或 NDVImax。TLNDVI 离 散 分 布 的 2 个 不 同 范围 分 
别 对 应 了 研究 区 内 的 高 山区 域 和 低地 平原 区 。 高 
海拔 地 区 气温 低 .冬季 长 .生长 季节 短 , 如 图 7bl 和 
图 7b2 中 的 A2 和 B2 区 ,因此 , 相 较 于 生长 季节 较 长 
的 区 域 , 即 使 其 NDVImax 可 能 相同 或 相近 ,TLNDVI 
却 要 小 的 多 (图 7b1)。 然 而 ,山区 具有 降水 丰富 的 
优势 ,夏季 温度 升 高 ,低温 限制 被 解除 后 , NDVImax 
则 远 远 高 于 干旱 的 低地 车 漠 区 (图 7b2)。 

NDVImax 和 TINDVI 均 呈现 明显 的 海拔 分 异 。 
水 和 热 是 控制 植被 空间 分 布 的 基本 因素 。 然 而 水 
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注 :Al 和 A2 区 域 的 散 点 NDVImax 相同 ,TI-NDVI 不 同 ,Al 的 TINDVI 大 于 A2 的 TINDVI( 红 色 和 矩 形 , 红 色 曲 线 );Bl1 和 B2 区 域 的 散 点 TI-NDVI 相 
同 ,NDVImax 不 同 ,B1 的 NDVImax 小 于 B2 的 NDVImax( 蓝 色 和 矩形 , 蓝 色 曲线 )。 
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图 7 NDVImax 和 TILNDYVI 的 散 点 图 (a) 及 NDVImax 或 TINDVI 相 同 但 TLNDVI 或 NDVImax 不 同 的 区 域 的 NDVI(b) 
Fig. 7 Scatter plot of NDVImax and TI-NDVI (a) and NDVI curves of pixels with same NDVImax but different TI-NDVI, NDVI 
curves of pixels with same TI-NDVI but different NDVImax (b) 


和 热 沿 海拔 梯度 变化 。 北 疆 地 人 处 亚 欧 大 陆 干 旱 区 
腹地 ,准噶尔 岔 地、 伊犁 河谷 和 塔 城 贫 地 等 海拔 较 
低 的 区 域 植被 受 缺 水 和 高 温 制约 ,所 以 这 些 区 域 的 
NDVImax 和 TI-NDVI 均 很 小 。 随 海拔 的 升 高 ,降水 
增加 ,气温 降低 ,对 植被 生长 的 限制 逐渐 减缓 ,ND- 
VImax 和 TI-NDVI 均 逐步 增 大 。 当 海拔 达到 一 定 高 
度 且 水 和 热 均 最 为 适宜 时 ,NDVImax 和 TI-NDVI 达 
到 最 大 ,如 天 山 和 阿尔 泰山 的 半山 腰 等 地 区 。 随 着 
海拔 的 升 高 ,植被 受 低温 和 饱和 土壤 水 的 限制 ,ND- 
VImax 和 TINDVI 又 开始 下 降 , 如 天 山 和 阿尔 泰山 
的 高 海拔 地 区 。 海 拔 对 水 热 条件 的 再 分 配 是 草地 
NDVImax 和 TI-NDVI 空 间 分 异 的 根本 原因 ,这 也 是 北 
性 草地 NDVImax 和 TILNDVI 均 呈 海 拔 分 异 的 原因 。 
4.2 NDVImax 与 TI-NDVI 变化 差异 的 解释 
NDVImax 和 TI-NDVI 在 年 际 波动 (图 3) 和 变化 
趋势 的 空间 格局 (图 4c) 存 在 差异 。 用 生长 季 内 不 
同时 间 的 NDVI 绘 制 的 曲线 不 仅 显 示 了 NDVImax 
(CNDVI 曲 线 的 顶点 ) ,而 且 还 显示 了 TI-NDVI( 遇 线 
上 累积 的 NDVI) 的 大 小 (图 8)。 因 此 ,通过 比较 不 


同 生 长 季 或 不 同 地 区 的 NDVI 曲 线 ,可 以 揭示 NDVL 
max 与 TI-NDVI 变 化 差异 的 原因 。 

对 于 整个 研究 区 域 ,2000 一 2001 4E , 2005— 
2006 年 、2014 一 2015 年 和 2018 一 2019 年 的 区 域 平 
35 NDVImax fI TI-NDVI 5t i8 5c 892 35 (F8 3)... 
图 8al 可 以 看 出 ,2001 年 NDVI 曲 线 的 峰值 较 低 ,但 
2001 年 第 209-305 d HY NDVI Lt 2000 年 大 得 多 , 表 
Hj 2001 年 仲夏 到 秋季 的 草地 状况 较 好 。 这 是 
2000—2001 4 NDV Imax 下 降 而 TI-NDVI 增 加 的 原 
因 。 从 图 8a2、 图 8a3 和 图 8a4 中 可 以 看 出 ,2005 年 、 
2015 年 和 2019 年 NDVI 曲 线 的 峰值 均 大 于 2004 年 、 
2014 年 和 2018 年 ,但 春秋 季 的 NDVIs 均 较 小 。 
2004—2005 年 , 2014—2015 年 和 2018—2019 年 ， 
NDVImax 增 加 ,TI-NDVI 减 少 。 

对 于 NDVImax Fil TI-NDVI 变化 趋势 的 空间 差 
异 。 天 山 山脉 和 塔 城 舍 地 周边 的 大 部 分 地 区 NDVI- 
max 降低 而 TI-NDVI 增 加 。 相 反 , 阿 尔 泰山 大 部 分 
地 区 NDVImax 增 加 而 TI-NDVI 降 低 ( 图 4c)。 对 比 
人 研究 期 前 5 a(2000—2005 年 ) 与 后 5 a(2015—2019 
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注 : 图 bl: NDVImax 减 小 .TINDVI 增 大 ;图 b2: NDVImax 增 大 、TI-NDVI 减 小 ;图 b3: NDVImax 增 大 、TI-NDVI 增 大 。 


图 8 TI-NDVI 5 NDVImax 变化 趋 热 相反 年 份 及 两 者 变化 趋势 不 同 铸 加 类 型 在 研究 共 


前 期 和 末期 NDVI 生 长 季 内 的 变化 


Fig. 8 Changes of TI-NDVI and NDVImax in the years with opposite trends and different superimposed types of the two trends in 


the NDVI growing season at the beginning and end of the study period 


年 ) 的 平均 曲线 (图 b) 可 知 ,在 天 山 及 塔 城 贫 地 周边 
山区 , 2000 一 2004 年 NDVI 曲 线 的 峰值 比 2015 一 
2019 年 略 大 ,但 2015 一 2019 年 曲线 中 第 65~145 d 的 
所 有 NDVI 均 大 于 2000 一 2004 年 (图 8pb1) 的 值 ,这 
表明 尽管 NDVI 峰 值 有 所 下 降 ,但 春季 NDVIs 则 有 
所 改善 ,提高 了 这 些 地 区 的 TI-NDVI。 对 于 阿尔 泰 
山地 区 ,夏季 NDVI 虽 有 增加 ,但 春季 和 秋季 减少 

多 (图 8b2) ,导致 NDVImax 增 加 而 TI-NDVI 减 少 。 
图 8b3 为 准噶尔 盆地 广阔 洼地 的 平均 曲线 ,其 中 
NDVImax fll TI-NDVI ZE 2000—2019 年 间 均 呈现 增 


加 趋势 (图 4c)。 可 以 看 出 ,整个 生长 季 的 NDVI 均 
有 所 提高 (图 8b3)。 

单一 生长 季 的 NDVI 曲 线 代表 了 草地 的 年 生长 
过 程 。 这 种 一 年 一 度 的 生长 过 程 通常 在 时 间 和 空 
间 上 都 有 所 不 同 。 而 年 生长 过 程 的 年 际 变化 也 是 
植被 动态 的 重要 特征 之 一 ,但 这 一 点 在 过 去 的 草地 
动态 监测 中 则 很 少 被 考虑 。 与 NDVImax 4H EE , TI- 
NDVI 积 累 了 草地 整个 生长 季 的 NDVI, 这 使 得 TI- 
NDVI 成 为 在 表征 草地 年 生长 过 程 方面 具有 价值 的 
量化 指标 。 
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5 结论 


(1) 北 疆 地 区 NDVImax 和 TI-NDVI 的 空间 变化 
基本 相同 , 均 表 现 出 明显 的 海拔 分 异 。TILNDVI 总 
体 随 NDVImax 的 增 大 而 增 大 ,但 NDVImax 或 TI-ND- 
VI 相同 的 区 域 ,其 TI-NDVI 或 NDVImax 却 存在 较 大 
差异 。 

(2) 2000—2019 年 北 疆 地 区 草地 TI-NDVI 和 
NDVImax 的 总 体 呈 现 增加 趋势 (P<0.01) ,但 空间 
上 ,整个 北 疆 17.55% 草 地 的 TI_NDVI 和 NDVImax 呈 
现 相反 的 变化 趋势 , 旦 多 位 于 草地 覆盖 度 较 高 的 山 
区 。 尤 其 是 阿尔 泰山 与 伊犁 河谷 山区 ,阿尔 泰山 草 
地 主要 表现 为 NDVImax 增 加 而 TLNDVI 减 小 ,伊犁 
河谷 山区 草地 则 主要 表现 为 NDVImax 减 小 而 TI- 
NDVI 增 加 。 

(3) 在 时 空 异 质 性 上 , 北 疆 阿 尔 泰 山 \ 伊 犁 河谷 
及 塔 城 贫 地 周边 山区 等 草地 高 覆盖 区 ,TI-NDVI 的 
Cys 和 Cv 均 高 于 NDVImax,TLNDVI 能 弥补 草地 动 
态 检 测 中 NDVI 交 饱和 缺陷 的 影响 ,能 更 准确 地 表 
达 高 覆盖 草地 的 时 间 异 质 性 。 
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Spatio-temporal heterogeneity evaluation of grassland TI-NDVI and 
NDVImax in northern Xinjiang 


JIAO Ayong, CHEN Fulong, YAN Junjie, LING Hongbo**, SHEN Ruihua' 
(1. College of Water Conservancy & Architectural Engineering, Shihezi University, Shihezi 832003, Xinjiang, China; 
2. Institute of Resources and Ecology, Yili Normal University, Yining 835000, Xinjiang, China; 
3. Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, Xinjiang, China; 
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Abstract: Grassland change is an important component of global change, which has attracted considerable 
attention. The temporal and spatial heterogeneity of grassland dynamics is the main concern in evaluating 
grassland dynamics. Northern Xinjiang, which is characterized with diverse grassland types, was selected as the 
research area. In this study, we calculated the time-integrated normalized vegetation index (TI-NDVI) and annual 
maximum NDVI (NDVImax) on the basis of the MODIS NDVI data. Using spatial analysis technology of GIS, 
mathematical statistical methods of coefficient of variation (Cy), and Mann-Kendall non-parametric statistics, the 
dynamic changes of grassland in northern Xinjiang were analyzed from 2000 to 2019, and the comparative 
advantages of TI- NDVI and NDVImax in expressing the temporal and spatial heterogeneity of grassland were 
explored. Results indicated that (1) the grasslands in northern Xinjiang, characterized by NDVImax and TI- 
NDVI, showed evident altitudinal differentiation. In general, the TI-NDVI increased with the increase of NDVImax. 
However, the areas with the same NDVImax or TI-NDVI showed great differences in TI-NDVI or NDVImax. (2) 
From 2000 to 2019, the grassland TI- NDVI and NDVImax in the northern Xinjiang showed a significant 
increasing trend (P<0.01), but the spatial differentiation of the changing trends of TI-NDVI was different from that 
of NDVImax. 17.55% of the grassland in northern Xinjiang showed opposite changing trends in TI- NDVI and 
NDVImax. For Altai Mountains and the mountains around Ili Valley, which are characterized with grassland of 
high coverage, the NDVImax and TI-NDVI showed opposite changing trends. (3) The Cy of TI-NDVI was higher 
than NDVImax in temporal and spatial dimensions in grassland with high coverage in northern Xinjiang. 
Furthermore, TI- NDVI was more sensitive to the temporal and spatial heterogeneity of high-coverage grassland, 
which can weaken the influence of saturation defect of NDVI in grassland dynamic evaluation to a certain extent. 


Keywords: spatio-temporal heterogeneity; TI-NDVI; NDVImax; grassland; northern Xinjiang 


